arXiv:1504.07349v3 [cond-mat.mtrl-sci] 18 May 2015 


Surface-termination dependent magnetism and strong perpendicular magnetocrystalline anisotropy 

of a FeRh (001) thin film: A density-functional study 


Soyoung Jekal*, S. H. Rhim'fl Soon Cheol Hong'J^ Won-joon Son^, and Alexander B. Shick^ 

* Department of Physics and Energy Harvest-Storage Research Center, 

University of Ulsan, Ulsan, 680-749, Republic of Korea 
^Samsung Advanced Institute of Technology, 

Suwon, 443-803, Republic of Korea 
^Institute of Physics, ASCR, Na Slovance 2, 

CZ-18221 Prague, Czech Republic 

(Dated: May 19, 2015) 

Magnetism of FeRh (001) films strongly depends on film thickness and surface terminations. While magnetic 
ground state of hulk FeRh is G-type antiferromagnetism, the Rh-terminated films exhibit ferromagnetism with 
strong perpendicular MCA whose energy +2.1 meV/D is two orders of magnitude greater than bulk 3d magnetic 
metals, where □ is area of two-dimensional unit cell. While Goodenough-Kanamori-Anderson rule on the 
superexchange interaction is crucial in determining the magnetic ground phases of FeRh bulk and thin films, 
the magnetic phases are results of interplay and competition between three mechanisms the superexchange 
interaction, the Zener-type direct-interaction, and magnetic energy gain. 

PACS numbers: 68.37.Ef,75.70.Tj,75.70.Rf 


FeRh alloys have attracted significantly because of their 
various intriguing physics phenomena includi^ magneto¬ 
caloric effect and huge magnetoresistancellUj]. Transi¬ 
tion between antiferromagnetism (AFM) and ferromagnetism 
(FM) occurs above room-temperature about 350 K and ultra¬ 
fast phase transition of magnetic phases in the FeRh alloys is 
induced by femto-second laser, which have drawn more atten¬ 
tion due to possibility of applications for heat-assisted mag¬ 
netic recording (HAMR)|l5l-l31- 

Furthermore, feasibility of room-temperature memories 
based on the AFM spintronics has been successfully demon¬ 
strated in the FeRh alloys very recently utilizing anisotropic 
magneto-resistance (AMR) of the bistable AFM stateslIsl-U^. 
This AFM^intronics has some advantages over that based on 
FM states 11814 1611 owing to the absence of stray magnetic field 
from the zero net magnetization and the insensitivity to the 
external magnetic fields. Regardin g Fe Rh thin films, several 
experiments have been rer)orted il7l4l9|] . Instead of the G-type 


antiferromagnetic (G-AFM) in bulk, some thin films exhibit 
FM states stabilized at the interface with metal, while the FM 
states are unstable at interface with oxide ifl^ . Interestingly, it 
was reported that an electric field of only a few volts is neces¬ 
sary to drive the AFM-FM transition for the epitax ially grown 
FeRh films on the ferroelectric BaTi03 substrate 1 1811 . Also, 
it was revealed that the spin orientation of the FeRh film on 
the MgO (001) depends on the strength of lattice strain and 
magnetic state 1 191] . 

In this paper, magnetism and magnetocrystalline anisotropy 
(MCA) of FeRh films are investigated using first-principles 
calculations. It is found that magnetism and MCA are sig¬ 
nificantly affected not only by film thickness but also by the 
surface-terminations. The Rh-terminated films are more sta¬ 
ble in FM state by quite big energy differences relative to 
the magnetic ground G-AFM state in bulk. Furthermore, 
the Rh-termination exhibits strong perpendicular MCA of 


+2.1 meV/D, where □ is two-dimensional (2D) unit cell 
area. The Fe-terminations, on the other hand, are in G- 
AFM states as in bulk. The strikingly different behavior of 
these two terminations is explained mainly in the framework 
of Goodenough-Kanamori-Anderson (GKA) rule on the su¬ 
perexchange interaction ll^ 211 with the Zener-type direct ex¬ 
change interaction i22ll also taken into account. 

Density functional calculations are performed using Vi¬ 
enna ab initio Simulation package (VASPl lEill . Results by 
VASP, particularly MCA, have been double-checked with the 
highly precise full-potential linearized augmented plane wave 
(FLAPW) method 12411 . Generalized gradient approximation 
(GGA) within projector augmented-wave (PAW) scheme ll^ 
is employed for the exchange-correlation interaction. k 
meshes of 24x24x24 and 24x24x1 in Monkhorst-Pack 
scheme are used for bulk and films, respectively. For wave 
function expansions, 500 eV is used for cutoff energy. Con¬ 
vergence with respect to cutoff energy and number of k points 
are seriously checked. Magnetic structures of the bulk FeRh 
are illustrated schematically in Fig. l(a)-(c) for A-, C- and 
G-type AFM, respectively, where arrows represent the direc¬ 
tion of magnetic moments of Fe atoms. To account for the 
AFM structures, a tetragonal magnetic unit cell, c(2 x 2) in 
the xy plane and doubled along the z axis, is taken. To in¬ 
vestigate thickness and surface-termination dependent mag¬ 
netism of the FeRh (001) films, two terminations, Fe- and Rh- 
termination with thickness from 3- to 15-ML have been taken 
into account, where we allow relaxation along the z-axis while 
fixing two-dimensional (2D) lattice constant as the calculated 
bulk value (3.007 A) of the G-AFM. 

First, magnetism of bulk FeRh is presented. AFM in G-type 
is found to be most stable from total energy calculations by 
48.3, 56.7, and 187 meV/Fe with respect to FM and AFM in 
A- and C-type, respectively. Since the C-AFM has the high¬ 
est energy, it will be excluded in forthcoming discussion if 
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FIG. 1: Schematic diagrams of magnetic structures of (a) A-, (b) C- 
and (c) G-AFM states of the bulk FeRh, and (d) exchange interac¬ 
tion between Fe atoms in the G-AFM state. Red and blue spheres 
at the corners represent Fe atoms with different magnetic states, re¬ 
spectively, while grey spheres represent Rh atom. In (d), cylinder 
connecting Fei and Fe 2 denote the 180° superexchange, while red 
cylinders connecting Fei-Rh-Fe 3 do the 90° superexchange. Blue 
dotted line which connects Fei-Fe 4 denotes the direct exchange in¬ 
teraction. 


not necessary. Calculated lattice constant (3.012 A) of the 
bulk FeRh in FM state is a little bit larger than that in the G- 
AFM (3.007 A). Magnetic moments of Fe and Rh atoms in 
the G-AFM (FM) are 3.158 (3.144) Hb and 0.00 (1.041) iig, 
respectively. Lattice constants and magnetic moments in this 
work are reasonably consistent with experiments ll^ Wh and 
previous calculation ll^ . 

In FeRh (001) films, the interlayer spacing exhibits oscil¬ 
latory feature; the topmost surface moves downward whereas 
the subsurface layers do upward, and so on, where the inner¬ 
most layer converge to bulk behavior [See Supplementary In¬ 
formation for fully relaxed interlayer spacings of the Fe- and 
the Rh-terminated film in Fig. SI (a) and (b), respectively]. 7- 
or 9-ML films are thick enough due to short metallic screening 
length. 

Total energy differences of the FM (A-AFM) and the G- 
AFM, which we denote as AE = Efm{Ea) — Eq, are pre¬ 
sented in Fig. 2(a) and (b) for the Fe- and Rh-termination, 
respectively, where FM, A, and G in subscripts stand for FM, 
A-AFM, and G-AFM, respectively. As shown in Fig. 2(a), 
the G-AFM is the most stable for the Fe-termination regard¬ 
less of thickness. We recall here that magnetism is so deli¬ 
cate physics phenomenon that it can behave differently in re¬ 
duced dimension ll^ . As film gets thinner, AE decreases 
while AE=36.6 meV/Ee for 9-ML is still less than that of bulk 
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FIG. 2: The thickness dependent total energy differences of the FM 
and A-AFM states from that of the G-AFM state, AE = Efm{Ea) — 
Eq ; (a) the Fe-terminated and (b) the Rh-terminated films. 


(48.3 meV/Ee). As shown in Eig. 2(b), magnetism of the Rh- 
termination is totally different from the bulk. EM is favored 
over G-AEM up to 15-ML by big energy difference whose 
absolute value is larger than 70 meV/Ee for 7-ML. Notewor¬ 
thy, there is a crossover from EM to G-AEM when the number 
of layers exceeds 17-ML. Eurthermore, even A-AEM is more 
stable than G-AEM for films thinner than 7-ML. Magnetic 
moments of Ee and Rh atoms in both terminations are listed in 
Table 1 and 2 for their respective magnetic ground states, G- 
AEM and EM. The surface layers have almost the same mag¬ 
netic moments as the bulk unlike other surfaces of magnetic 
elements, Ee, Co, and Ni, where moments are strongly en¬ 
hanced. 

Before we discuss magnetism the EeRh films, we provide 
here a detailed analysis of bulk magnetism. Cg and t 2 g are 
degenerate in G-AEM and EM due to the cubic symmetry, 
whereas those in A- and C-AEM are no longer degenerate 
owing to the tetragonal magnetic structure (see SI Eig. S2). 
Strong hybridization between Ee and Rh d states brings in 
almost fully occupied (unoccupied) majority (minority) spin 
states of Ee d orbital, particularly Cg states for all magnetic 
phases, which results in enhanced moments of Ee atoms. The 
nearly half-filled band favors AEM as in bulk Mn and Cr. It 
is noteworthy that the majority spin bands of Rh are almost 
fully occupied in the EM similarly to Ee, but featureless in 
AEM states. 

The magnetic phase of EeRh alloy can be explained 
by interplay between three mechanisms - superex¬ 
change interaction i^ 21], Zener-type direct-exchange 
interaction!^], and magnetic energy gain. In the frame¬ 
work of GKA rule, whether EM or AEM is preferred by 
superexchange interaction is explained by magnetic ion- 
ligand-magnetic ion angle. Here we view Ee atoms as 
magnetic ions and delocalized s and p orbitals of Rh as ligand 
orbitals in GKA rule. Eig. 1(d) schematically illustrates 
magnetic interactions between Ee atoms to be involved in 
determining magnetic structure. The GKA superexchange 
interactions are shown as solid lines and Zener-type direct 
interaction as a dotted line. In accord with GKA rule on 
magnetic coupling, Eei prefers AEM coupling to Ee 2 and 
EM coupling to Ee 3 because angles of Eei-Rh-Ee 2 and 
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FIG. 3: MCA energies, defined as Umca = ^(~^) of the 

FeRh (001) thin films for (a) Fe- and (b) Rh-terminated films, respec¬ 
tively. Solid (dashed) lines denote MCA of FM (G-AFM) phase. 


Fei-Rh-Fe 3 are 180° and 109.5° (which is close to 90°), 
respectively. 

On the other hand, the interaction between Fei and Fe 4 is 
more or less direct since Rh atom is not much involved in this 
coupling. The half-filled eg states directing along the principal 
axes are more involved in the Zener direct interaction between 
Fei and Fe 4 compared to the t^g states, which results in AFM 
coupling j^ . Since the d states are highly localized giving 
little wave function overlap, the Zener-type direct interaction 
between Fei and Fe 4 must be weak. As a result of combina¬ 
tion of the superexchange and the direct interactions discussed 
above, G-AFM is most stable among other AFM states. In the 
FM states, on the other hand, there is magnetic energy gain 
due to the considerable magnetic moment of Rh atom (1.041 
/Tb), which reduces total energy to a certain degree, hence a 
FM state is more stable than A-AFM and C-AFM states even 
though it is less stable compared to the ground G-AFM state. 

In the Rh-termination, the 180° superexchange interaction 
disappears because of the absence of Fe layer above the Rh- 
terminated surface. The 90° superexchange interaction makes 
the A-AFM more stable than G-AFM. Instead, the magnetic 
energy gain of the surface Rh atom plays a key role in sta¬ 
bilizing FM. Hence, FM is the magnetic ground state in the 
Rh-terminated FeRh thin films. In the Fe-termination, on the 
other hand, G-AFM is most stable as in bulk since the 180° 
superexchange interaction still works. 

In Fig. 3, calculated MCA energies, Emca = E (—>^) — ^(t), 
are presented as function of thickness for the Fe- and the Rh- 
terminated films in their for G-AFM and FM states. Emca^ 
for other magnetic states are also shown for comparison. 
From the definition of E^ca, positive (negative) value im¬ 
plies perpendicular (parallel) magnetization to the surface nor¬ 
mal. Interestingly, the Rh-termination in FM state shows quite 
strong persistent perpendicular MCA regardless of thickness, 
whereas the Fe-termination exhibits in-plane MCA for all 
magnetic states. In particular, Emca =+2.1 meV/D of the Rh- 
termination is greater than 3d magnetic metals such as bulk 
Fe, Co, and Ni by two orders of magnitude and larger than 
magnetic films by several factors. From the fact that both the 
Fe-terminated FM states and the Rh-terminated G-AFM state 
show parallel MCA, the magnetic states are not a key factor in 
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FIG. 4: Projected density of states of d orbitals of 9-ML FeRh film 
in the Fe-termination (upper panels) and the Rh-termination (lower 
panels) for FM and G-AFM states. 5 and 5-1 indicate the top surface 
and subsurface layer, respectively. 


determining MCA, but the surface termination is. Moreover, 
the saturated feature of Emca =+2.1 meV/D with respect to the 
thickness implies that magnetism of the Rh surface plays a key 
role in determining the strong perpendicular MCA. To discuss 
the role of the thickness and the surface termination on mag¬ 
netism as well as MCA, we present density of states (DOS) 
of the Fe- and the Rh-terminated 9-ML film is presented in 
Fig. 4, where left and right columns represent DOS’s of Fe 
and Rh atoms for better comparison as in Fig. S2. First and 
second rows are for the Fe-terminated film in G-AFM and FM 
states; third and forth rows are for the Rh-terminated film in 
FM and G-AFM states. DOS features do not differ very much 
from the bulk, despite lifted degeneracies in eg and tig states 
at the surfaces and some changes in surface atoms. In particu¬ 
lar, DOS’s of subsurface atoms are essentially the same as the 
bulk. This result confirms that the surface atoms indeed play 
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the key role in determining magnetism in film geometry, as 
saturated Emca with respect to thickness implies [see Fig.3]. 

To confirm and analyze calculated results on MCA, ad¬ 
ditional calculations have been carried out using FLAPW 
method using GGA for 5-, 7-, and 9-ML Rh-terminated films. 
Results by FLAPW method are quite consistent with VASP 
method; i) FM is much more stable than G-AFM states, 
ii) magnetic moments of the surface Rh atoms are 1.107, 
1.106, and 1.107 (Compare with SI Table S2) for 5-, 7-, and 
9-ML, respectively, and iii) Emca =+L97, h- 2.41, and 2.10 
meV/D (cf. Fig. 3), respectively. Emca^ are further de¬ 
composed into individual atomic contribution and different 
spin-channels lIiHl . It is found that the surface Rh atom with 
stronger spin-orbit cour)Iing |l33 - [^ and the fj, channel play 
a dominant role in determining MCA. Results using FLAPW 
adopting conventional von Barth-Hedin local spin density ap¬ 
proximation (LSDA) are also listed in SI Table S3. 

The role of the surface Rh layer in determining perpendic¬ 
ular MCA is well manifested in DOS: the lifted degeneracy of 
Cg states makes peak from more prominent in the minority 
spin band, which is just above the Fermi level. This peaked 
state contributes significantly to perpendicular MCA through 
(yz/zx,t \Ex\^,\) matrix i3ill in the Rh-termination. 

In order to elucidate the magnetic behavior when the 
FeRh(OOl) films are under an applied magnetic field or heat, 
we present total energy of the Fe-terminated FeRh(OOl) thin 
film as a function of angle between magnetic orientations of 
the closest Fe atoms.[See SI Fig. S3]. The energy barrier of 
the AFM-to-FM transition under an applied magnetic field is 
reduced with decreasing of thickness. This information might 
be useful in designing spintronics devices such as an AFM 
memor ylUlIi and HAMrII-H]. 

In summary, magnetism of the Fe- and the Rh-terminated 
FeRh (001) are studied for various film thicknesses. The ori¬ 
gin of stability of G-AFM in bulk is well explained in the 
framework of Goodenough-Kanamori-Anderson rules on the 
super-exchange interaction, where subsidiary Zener-type di¬ 
rect exchange interaction and magnetization energy are also 
taken into account. The thickness and the surface termina¬ 
tion turn out really significant as the two terminations give 
different magnetic ground state. The Fe-termination is sta¬ 
bilized in G-AFM as in bulk, while the energy difference 
between G-AFM and FM is greatly reduced. On the other 
hand, the Rh-termination strongly prefers FM when films are 
thinner than 15-ML. Furthermore, the Rh-termination in FM 
state exhibits quite strong perpendicular MCA, 2.10 meV/D 
, which is greater than bulk 3^f convention magnetic metals 
by two orders of magnitude and larger than 3^f magnetic films 
by several factors. Utilizing the thickness and the surface- 
termination dependent magnetism of FeRh film, highly desir¬ 
able materials can be tailored as required by spintronics de¬ 
vices. 
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